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Abstract 

Energy band structures are calculated for the new superconductor MgB2 and 
the related compounds by using the LDA and an FLAPW method. It is found 
that the strong three dimensional network in low-lying vr bands brings about two 
dimensional a holes in MgB2, which should be responsible for the superconductivity. 
The de Haas- van Alphen frequencies and the cyclotron masses are obtained for the 
Fermi surfaces. The a hole is not found in LiBC due to the less three dimensional 
network. MgB2C2 is found a semiconductor, but the top of the valence bands are 
similar to the a bands in MgB2 . The total energy calculations are also performed for 
Mgi_2;Al2:B2 to investigate the structural phase instability, experimentally observed 
around x =10%. 
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1 Introduction 



Prom the discovery of the superconductivity in MgB2[l], a lot of theoretical 
and experimental researches have been performed in order to reveal the mecha- 
nism of the really high superconducting transition temperature (Tgc) in metals 
and find a route to a new superconductor with much higher Tgc- It has been 
indicated that MgB2 is a phonon-mediated BCS s-wave superconductor by 
theoretical calculations [2, 3], the NMR study[4] and the high-resolution pho- 
toemission study [5], then detailed phonon dispersion is calculated [6, 7]. On the 
other hand, multiple superconducting gaps are also indicated by the specific 
heat measurement [8] , the tunneling spectroscopy study [9, 10] and the high- 
resolution photoemission study[ll]. In the electron doped system Mg1_a.Ala.B2, 
the structural instabilityhas been observed around x 10%, resulting in the loss 
of superconductivity. [12] 

In this paper, the Fermi surface property in MgB2 is theoretically obtained for 
the further microscopic measurements of single crystals. The bandstructures 
for iso-electronic compounds LiBC and MgB2C2 are calculated for the com- 
parison. Total energy calculations for Mgi_a:Ala:B2 are performed to clarify the 
structural instability. 



2 Method of Calculation 

Band structure calculations are carried out by using a full potential APW 
(FLAPW) method with the local density approximation (LDA) for the ex- 
change correlation potential. For the LDA, the formula proposed by Gunnars- 
son and Lundqvist[13] is used. For the band structure calculation, we used the 
program codes; TSPACE[14] and KANSAI-99. 

We used in the calculations the experimentally observed lattice constants: 
a= 3.083A and c=3.527A for MgBa. Muffin-tin (MT) radu are set as 0.4979 
a and 0.2742 a for Mg and B, respectively. Core electrons (Ne-core for Mg, 
He-core for B) are calculated inside the MT spheres in each self-consistent 
step. The LAPW basis functions are truncated at |k + Gj| < 3.33 x 27r/a, 
corresponding to 225 LAPW functions at the P point. The sampling points 
are uniformly distributed in the irreducible 1/24 th of the Brillouin zone, 333 
k-points (divided by 12, 12, 12) for both the potential convergence and the 
final band structure. The similar parameters are used in calculations for LiBC 
and Mgi_^Al^B2. 

MgB2C2 crystallizes in a base centered orthorhombic structure containing 8 
formula units in the primitive unit cell. In this case, calculations are performed 
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by using about 1,600 LAPW basis functions at the uniformly distributed 39 
k-points for the potential convergence and 369 k-points (divided by 8, 8, 8) 
for the final band structure in the irreducible 1/8 th of the Brillouin zone. 



3 Bandstructure and Fermi surface of MgB2 

The calculated energy bandstructure and the density of states are shown in 
Fig. 1 and Fig. 2, respectively. They are very similar to the previous calcula- 
tions[2,3], though Fig. 1 shows bandstructures along all the symmetry axes; 
the F-K-M-F lines is parallel to the A-H-L-A lines on the basal plane and the 
A, U and P axes are along the c-axis (see Fig. 3 (a)). 

Let us see the dispersion along the P axis; the dispersive band crossing the 
Fermi level. The band is originated from B-pz{7r) band and doubly degen- 
erated. These B-p^(7r) states arc non-bonding along the P axis within the 
honeycomb lattice, and anti-bonding (bonding) at the K (H) point for the 
inter-layer (along the c-axis) mixing. While B-p^y{a) states arc split to bond- 
ing and anti-bonding sates within the honeycomb lattice, then the top of the 
bonding states are located just above the Fermi level along the A axis. While 
the bottom of the anti-bonding bands are located around 1.1 Ry. The energy 
and its components foe the states near the Fermi level are listed in Table 1. 
Eventually, there arc four Fermi surfaces; two cylindrical hole Fermi surfaces 
along the A axis originating in the a band, and from the tt bands one hole 
around the T' axis one electron along S' axis, as shown in Fig. 3. The total 
number of holes and electrons are compensated. 

The unique characteristics in the bandstructure of MgB2 is the existence of 
the cylindrical cr-holes (Fig. 3 (a) and (b)). The situation has been discussed 
as the relative shift of the a and tt bands[2]. In fact, the totally non-bonding a 
states (about 0.91 Ry) is much higher in energy than the non-bonding vr state 
(0.65 Ry; the mid-point on the P axis). It should be ascribed to the ionic Mg^"*" 
strongly affecting the crystalline field for the B-p states. Moreover, here it is 
emphasized that the 3D network in the tt band is suggested by the dispersive 
band along the P axis. 

As was already discussed[2], graphite (primitive graphite C2) shows a 2D band- 
structure. In graphite, the inter-layer distance {diayer =3.348A) is smaller than 
in MgB2 {diayer =3.528A), but the planar bond length (dbond =1.418 A) is much 
smaller than in MgB2 {dbond =1-780 A). The planar bond length determines 
the total band width in each case, because of the much smaller distance and 
more coordination in plane, though it depends on how extended the wave 
functions are. As a result, the valence band width is narrower in MgB2 than 
graphite. However, what should be noticed here is the ratio of diayer and dbond, 
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which is 2.36 in graphite and 1.98 in MgB2. The much smaller ratio in MgB2 
causes more 3D feature, resulting in the dispersive it band along the P axis, 
which obtains both the hole and electron Fermi surfaces. The a bands keep 
the 2D character, because the p^y orbitals are extended in a layer. Eventually, 
the strong 3D feature in low-lying tt bands brings about 2D Fermi surface of 
the a band in MgB2, as shown in Fig. 3. The 3D feature and the low- lying in 
TT bands are ascribed to the existence of Mg^+ ion. 

The calculated specific heat coefficient 'jband is 1.73 mJ/mol K^, while the 
experimental value ^exp is 2.6 mJ/mol K^.[8] Therefore, the averaged electron- 
phonon mass enhancement factor A is derived as 0.50 from ^yexp — {i + ^)^band- 

The calculated angular dependence of the de Haas-van Alphen (dHvA) fre- 
quencies are shown in Fig. 4. The dHvA frequencies and the cyclotron masses 
are summarized in Table 2. Several theoretical studies have pointed out the 
importance of electron-phonon coupling in the a band for the superconduc- 
tivity, then the superconducting gap in MgB2 can be expected to be larger 
on the a Fermi surfaces than on the tt Fermi surfaces. It is also expected to 
be larger mass enhancement in the cyclotron masses in the a Fermi surfaces, 
though they have not yet observed so far, unfortunately. Recently MgB2 single 
crystals are grown by a several groups, [15-17] then the dHvA oscillations is 
expected to be detected soon. The experimentally observed cyclotron masses 
will be the direct evidence for the two superconducting gaps in MgB2. 



4 Bandstructures of LiBC and MgB2C2 

It is worth to see bandstructures of the related compounds. LiBC and MgB2C2 
are iso-electric compounds with MgB2, though the crystal structure is slightly 
different. [18] LiBC can be obtained, if Mg and one of B in MgB2 are replaced 
by Li and C. The stacking along the c-direction has not been experimentally 
determined, so we assume B-B-B...(C-C-C...) stacking. MgB2C2 crystallizes in 
a orthorhombic structure, as shown in Fig. 5, in which there are a distorted B- 
C honeycomb lattice layer and a Mg layer. In LiBC, diayer =3.53 A, dbond =1-59 
A and diayer/dbond =2.22. While, in MgB2C2, diayer =3.74 A, dbond =1-58 A 
in average and diayer /dbond =2.37. Therefore, both compounds are expected to 
be less 3D than MgB2. 

Figure 6 shows the bandstructure for LiBC. The space group of LiBC is P6m2, 

then the degeneracy along the P-axis is lifted. Li^+ ion has less attractive for 
the TT bands, then the bonding a bands get under the Fermi level. Eventually, 
there are hole and electron Fermi surfaces both from the tt bands, as shown 
in Fig. 7. 
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More interesting bandstructure is found in MgB2C2. As shown in Fig. 8, 
MgB2C2 becomes a semiconductor in consistent with the measured conductiv- 
ity. [18] The opening the gap may be due to the less anisotropic hybridyzation 
in the lower symmetry. However, the dispersion of top of the valence bands 
are very fiat along the A axis, which is similar to the a bands in MgB2. Cylin- 
drical Fermi surface, which is relevant in the superconductivity in MgB2, can 
be obtained with hole doping. 



5 Lattice Constants and Structural Instability in Mgi_a;Ala;B2 

With substituting Al for Mg in MgB2, the superconducting transition temper- 
ature gradually decreases and structural instability appears with Al content 
from 10 % to 20%. [12,19] After the instability, the superconducting transition 
temperature rapidly decreases. Partial Al substitution for Mg means electron 
doping to the MgB2 bandstructure, leading to vanishing a hole Fermi surfaces. 
One can expect the structural phase instability is related to the disappearance 
of the a hole. 

To reveal the structural instability, the total energy calculations with a several 
concentration and the lattice constant are performed. In the calculations, the 

lattice constant a is assumed to be constant, as its small Al concentration 
dependence (a=3.083A in MgB2 to 3.006 A in AIB2) and the virtual lattice 
approximation is adopted, in which the atom with the atomic number of 12. a; 
is used for the case of Al concentration x. 

The total energies are obtained at the points in Fig. 9, then the optimized 
lattice constant c for each concentration x are determined by using parabolic 
fitting. The total energy as a function of c is well reproduced as a parabolic 
function, so that the present calculations cannot predict the observed struc- 
tural instability for any Al concentration. There is another structural insta- 
bility reported around x =70%, where the superconductivity disappears. [20] 
Nevertheless the optimized c decreases monotonically with increasing x, as in 
Fig. 9. Concerning that Aluminum ordering is observed for x — 0.5[20], the 
virtual lattice approximation may not be applied to the system Mgi_3;Al3;B2. 

For AIB2, with using the measured a=3.006A, the optimized c is success- 
fully obtained as 3.245A for the experimental value 3.25lA. With the value 
a=3.083A for MgBa, which is used in Fig. 9, c=3.24lA for AIB2. In each case, 
the error is less than about 0.3% for AIB2. However, for MgB2, the optimized 
c is 3.457A which is 2% smaller than the measured value 3.527A. The lattice 
constants decrease in low temperatures, but c is reported as 3.51504(4)A at 
15K, which is only 0.3 % smaller than in room temperature. [21] This 2% error 
should make some sense and corresponds to the jumps of two discontinuities 
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(structural instabilities) around 20% and 70% Al concentration. The larger 
inter-layer distance (c) could decrease the number of the a holes. Therefore 
there may exist some mechanism beyond the LDA that favours the small 
number of the a holes and is related to the superconductivity. This should be 
solved in the future. 
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Table 1 

The energy and the component (%) inside the MT spheres for the states near the 
Fermi energy 0.6513Ry. They consist of mainly B-p component, but arc extended 
to Mg-site then Mg-p and -d components appear. "-" means the component is 
forbidden by the symmetry. 
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Table 2 



The calculated dHvA Frequencies F in unit of Tesla and the cyclotron masses in 
unit of free electron mass mo- 
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MgB, 




Fig. 1. The calculated bandstructure for MgB2. 
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Fig. 2. The calculated total and partial density of states for MgB2. 
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Fig. 3. The calculated Fermi surfaces for MgB2; (a) hole Fermi surface from the 3rd 
band, (b) hole Fermi surface from the 4th band and (c) electron Fermi surface from 
the 5th band. 12 
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Fig. 4. The calculated angular dependence od the dHvA frequencies for Mg 
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5. The crystal structure of MgB2C2; (a) side view and (b) top view. 
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LiBC 




r r K'ta = r*A 8 hB'L s amulk'h 

Fig. 6. The calculated bandstructure for LiBC. 
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5th band electon 



Fig. 7. The calculated Fermi surfaces for LiBC; (a) hole Fermi surface from the 4th 
band and (b) electron Fermi surface from the 5th band. 
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Fig. 8. The calculated bandstructure for MgB2C2. 
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Fig. 9. The optimized lattice constant c (open curcle) vs Al-concentration x, keeping 
a constant. The total energies are calculated at the point denoted by x. 
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